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' Abstract 

. 

Considering gravitational collapse of Type I matter fields, we prove that, given 
an arbitrary C 2 - mass function M(r, v) and a C 1 - function h(r, v) (through the cor- 
responding C 1 - metric function u(t,r)), there exist infinitely many choices of energy 
distribution function b(r) such that the 'true' initial data (M, h(r,v)) leads the col- 
lapse to the formation of naked singularity. We further prove that the occurrence of 
such a naked singularity is stable with respect to small changes in the initial data. 
We remark that though the initial data leading to both black hole and naked singu- 
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bJT). larity form a "big" subset of the true initial data set, their occurrence is not generic. 

The terms 'stability' and 'genericity' are appropriately defined following the theory of 
dynamical systems. The particular case of radial pressure p r (r) has been illustrated 
in details to get clear picture of how naked singularity is formed and how, it is stable 



with respect to initial data. 
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1 Introduction 

In recent papers [1, 2, 3] Joshi, Dwivedi and Goswami have discussed the role of initial data 
in spherically symmetric gravitational collapse for Type I matter fields. It is shown that 
given the density (p) and pressure profiles (p r ,Pe) a t the initial spacelike hypersurface from 
which the collapse evolves, there is a wide choice available to choose the velocity function 
and rest of the initial data so that the end state of collapse is either a black hole (BH) 
or a naked singularity (NS). This result is significant for two reasons:- (1) it produces a 
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substantially "big" initial data set which under gravitational collapse results into a naked 
singularity. (2) Type I matter fields include most of the known physical forms of matter 
like dust, perfect fluid. At least in spherically symmetric case, this result shows that while 
considering black hole physics, one must discard certain part of initial data set that leads 
the collapse to a naked singularity. However, all of the initial data is not independent, so 
it poses a question, which initial data set shall evolve into a NS as a result of gravitational 
collapse? Further, having a NS for certain initial data set is not enough because if it 
looses its characteristic for small perturbation in the neighbourhood of that certain initial 
data set then the NS is not serious enough to challenge cosmic censorship conjecture [4]. 
Therefore, the all important question ie. is the NS developed from certain initial data set 
a stable one? 

This motivates us to study both the above mentioned questions. We prove that given 
a C 2 -mass function M(t, r) and C 1 - velocity function v(t, r) on an initial spacelike hyper- 
surface, there always exists a C 1 - function b(r) such that the total initial data comprising 
of mass function (includes density and pressures), a function v{r) and energy distribution 
function b(r) evolves into a naked singularity as a result of gravitational collapse. 

F.I. Cooperstock and et al. [5] have shown that if the radial pressure p r is positive , 
r > and having strictly positive mass function then non-central singularity is covered by 
the horizon irrespective of the sign of tangential pressure pg. The stability of central shell- 
focusing singularity at r = that we are discussing in this paper is with respect to small 
changes in the initial data comprising of mass function and velocity function. "Small" 
changes in the initial data are with respect to appropriate norm on the the functions 
space of all physically reasonable initial data. This concept of stability is analogous to 
that of structural stability in the theory of dynamical systems and is explained in section 
4. Essentially, we prove that the set of initial data leading the collapse to central shell 
focusing naked singularity is an open subset of the set of all initial data in a suitable 
functions space. That we get a covered singularity which is non-central when p r > does 
not disturb the stability of central shell focusing singularity (naked or covered) with respect 
to small changes in initial data. 

In section 2, we briefly summarize analysis given in [3] and state the conditions on the 
initial data under which the collapse will lead to a naked singularity. In section 3, we 
analyze the condition of existence of naked singularity to prove our assertion mentioned 
above by using existence theory of first order ordinary differential equations. Thereafter, 
in section 4, we define and explain the concept of stability and genericity that we follow in 
this paper and we prove openness of initial data set. We, then, discuss genericity of this 
set. In section 5, we study the particular case p r (r) to demonstrate the involvement of 
p r and po as an initial data in the occurrence of NS, its stability and non-generic nature. 
Finally, we make some concluding remarks. 

2 Field equations, initial data and naked singularity 

In this section, we summarize analysis from [3]. The general spherically symmetric metric 
describing space-time geometry within the collapsing cloud can be described in comoving 
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coordinates (t, r, 9, 0) by 

dg 2 = _ e 2v(t,r) dt 2 + e 2Mt,r) dr 2 + R 2^ ^2 (1) 

where dfl 2 = d6 2 + sin 2 6 d(f> 2 is the metric on an 2-sphere. The stress-energy tensor for 
Type I field in a diagonal form is given by [6] 

T t t = -p,T r r =p r ,I* e = T+ = p e (2) 

The quantities p, p r and pe are the energy density, radial and tangential pressures respec- 
tively. We take the matter field to satisfy weak energy condition i.e. the energy density 
measured by any local observer be non-negative, so for any vector V 1 , we must have, 
T ik V' l V k > which means p>0; p + p r >0; p + p e >0. 
Einstein field equations for the metric (1) are 

F 

Pr = (4) 

v R 2 R V ; 

v , = 2(Pe ~ Pr) R' _ p' r ^ 

p + p r R P+Pr 

-2R> + R'^ + R^ = (6) 

G-H=l-^ (7) 

where we have defined G(t,r) = e~ 2 ^(R') 2 and H(t,r) = e~ 2u (R) 2 . 

The total mass in a shell of collapsing cloud of comoving radius r is represented by the 
arbitrary function F = F(t,r). The weak energy conditions imply F > 0. The regularity 
at the initial epoch t — U is preserved by F(ti, 0) = i.e. the mass function should vanish 
at the centre of the cloud. The density singularity in the space-time is at R = and R' = 0, 
the later one is due to shell- crossings and can possibly be removed from the space-time [7], 
so we consider here only a physical singularity where all matter shells collapse to a zero 
physical radius known as shell-focusing singularity. We can use scaling freedom available 
for the radial co-ordinate r to introduce the function v(t,r) by the relation 

R(t,r) = rv(t,r), (8) 

we have v(ti,r) = 1; v(t s (r),r) = and for collapse v < 0. Above relation is obtained by 
defining v(t,r) = R/r [8]. The time t = t s (r) corresponds to the shell-focusing singularity 
R — 0. The six arbitrary functions of the shell radius r as given by v(ti,r) = v {t), 
il>(U,r) = ip (r), R(ti,r) = r, p(U,r) = p (r), p r (t i: r) = p ro (r), p e (U, r ) = Pe a (r) evolve 
the dynamics of the initial data prescribed at the initial epoch t — ^. From equation (5), 
we obtain 

„. (r)= rf%z^__^_U (9 ) 

J0 \r{p + Pr ) Po + PrJ 
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which indicates that not all the above initial data are mutually independent and that z/ D (r) 
has the form, 

v {r) = r 2 g(r) (10) 

where g(r) is at least a C 1 function of r for r = 0, and at least a C 2 function for r > 0. 
Let us now assume the physically reasonable behaviour of the initial data at the center 
r = i.e. p f ro (0) = p'g o (0) = and p ro (0) — pe o (0) = 0. For details, we refer [1]. We 
have a total of five equations with seven unknowns, namely p,p r ,pg, u, ip, R and F giving 
us the freedom of choice of two functions. Selection of these two free functions, subject 
to the given initial data and the weak energy condition above, determines the matter 
distribution and the metric of the space-time and thus, leads to a particular time evolution 
of the initial data. The existence and uniqueness of solution of the system of field equations 
with above mentioned initial data has been discussed by Joshi and Dwivedi [9, 1]. The 
solution continues to exist in the neighbourhood of the singularity given by R = 0. A 
general mass function for the cloud can be considered as 

F(t,r) =r 3 M(r,v) (11) 

where M(r,v) is regular and continuously twice differentiable. Using equation (11) in 
equations (3) and (4), we obtain 

3M+r[M <r + M tV v'] M v 

P = TT~a A ; P r = 2 12 

v z {v + rv ) v z 

Then as v — > , p — > oo and p r — * oo i.e. both the density and radial pressure blow 
up at the singularity. The regular density distribution at the initial epoch is given by 
p {r) = 3M(r, 1) + rM(r, l) r . 

We take the initial surface to be the cloud given by < r < r c for some finite r c on 
which the initial data, namely the mass function F(t,r), the metric function v(t,r) and 
the function b(r) ( to follow) evolve as the collapse begins according to the Type I field 
equations. So, we take M(r,v) > and at least C 2 in V = [0, r c ] x [0,1], for varying 
r G [0, r c ], v E [0, 1]. Consider the general metric function, 

u(t,r) = A(t,R) (13) 

where A(t, R) is arbitrary, continuously differentiable function of t and R, with the initial 
constraint A(t i: R) = u (r). We use equation (13) in equation (5) to get 

Rv' 

2 Pe = RA, R (p + p r ) + 2p r + ^ (14) 

From above equation, we conclude that the tangential pressure also blows up at the sin- 
gularity. Also, the use of equation (13) in equation (6), yields 

G(t,r) = d(r-)e 2i - A -$ A ^ (15) 

where d(r) is another arbitrary continuously differentiable function of r. From equation 
(10), we generalize the form of A(t, R) as A(t, R) = r 2 gi(r, v) , where gi(r, v) is a continu- 
ously differentiable function and gi(r, 1) = g(r). Similarly, we have A — J A yt dt = r 2 g2(r,v) 
and at the initial epoch g 2 (r, 1) = g{r) . We write 

d(r) = l + r 2 b(r). (16) 
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where b(r) is the energy distribution function for the collapsing shells. Then using equa- 
tions (11), (13) and (15) in equation (7), we get 



y/RR = - e r 9l ^ v) 01 + r 2 b(r)]Re r2 ^(r,v) - R + r 3 M (17) 
where negative sign is chosen since, for the collapse, R < 0. Defining a function h(r,v) as 

e r 2 g 2 (r,v) _ i 

h(r,v) = = = 2g 2 (r,v) + 0(r 2 v 2 ) (18) 



r 



and substituting equation (18) in equation (17), we get 



yfiv = -^Je 2r2( -9i+g2) v b(r) + e 2r2 ^{vh + M). (19) 
Integrating the above equation , we have 

t(v,r)= V = (20) 

J v ^ e 2r 2 (9i+92) v b( r ) + e 2r2 si (vh + M) 

Here, the variable r is treated as a constant in the above equation. Expanding t(v,r) 
around the centre, we get 

t(v,r) = t(v,0) + rx(v) + 0(r 2 ) (21) 

where the function 

If' ^[vb'(0)+vh r (0,v) + MA0,v)} , , . 

X{ ] 2J V [vb(0)+vh(0,v) + M(0,v)} 3 / 2 ' 1 ' 

Now, it is clearly seen that the value of x(0) depends on the functions 6(0), M(0, v) 
and /i(0, v), which in turn, depend on the initial data at the initial surface t = ti. Thus, 
a tangent to the singularity curve t = t So is completely determined by the given set of 
density, pressure, velocity function v and function b(r). Further, from equation (19), we 
can write 



vv ' = x(v)yjvb(0) + vh(0, v) + M(0, v) + 0(r 2 ) (23) 

If the neighbourhood of the centre R = 0, r = gets trapped earlier than the formation 
of singularity, then it is covered ( i.e. occurrence of a black hole), and if otherwise happens 
the singularity is naked (ie. singularity can be observed locally or by a faraway observer) 
[7]. For examination of the nature of central singularity at R = 0, r = 0, we consider the 
equation for outgoing radial null geodesic, dt/dr = e^~ u ^ and test, if there would be any 
families of null geodesies coming out of the singularity. Further, we write the null geodesic 
equation in terms of the variables ( u = r a ,R), choosing (a = 5/3 ), and using equation 
(7), we obtain 



dR 




F 
R 



du 5 U JI \^/G[VG + y/H] 



(24) 



If the null geodesies terminate at the singularity in the past with a definite tangent, then 
at the singularity, we have dR/du > 0, in the (u,R) plane with a finite value. Hence, all 
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points r > on the singularity curve are covered since F/R — > oo with dR/du — > — oo and 
only the singularity at the centre r = could be naked. For the case, when R' > near 
the central singularity, we have 



, , R dR 
x Q = hm hm — = — — 

t-»t s r^o ii d u 



(25) 



and use of equations (23) and (24) yield x Q 3/2 = (3/2)^(0, 0)x(0) and the radial null 
geodesic emerging from the singularity in u) co-ordinates is R = x Q u, or in (t, r) plane, 
it is given by t — t s (0) = x Q r 5//3 . Therefore, if x(0) > then x Q > 0, and we get radially 
outgoing null geodesies coming out from the singularity , giving rise to a naked central 
singularity. However, if x(0) < then we have a black hole solution, as there will be no 
such trajectories coming out. For %(0) = , we will have to take into account the next 
higher order non-zero term in the singularity curve equation, and a similar analysis can be 
carried out by choosing a different value of a. 

Thus, it is clearly seen as to how the initial data determines the NS/BH phases as end 
states of collapse, since x(0) is determined by these initial profiles for the collapsing matter 
given by equation (22). 

3 Existence of energy distribution function b(r) lead- 
ing to NS 

In this section, we prove the assertion mentioned in the introduction. We choose b(r) to 
satisfy the differential equation on a constant v-surface 

l yfi[vV(r)+vh r (r,v) + M tr (r,v)] = 

2 [vb(r)+vh(r,v) + M{r,v)fl 2 K,) { 1 

for < r < r c where B(r,v) is a continuous function defined on V such that 

mn v] = l yfiW(0)v + vhAO,v) + M tr (0,v)] 
K,) 2 [vb(0)+vh(0,v) + M(0,v)} 3 / 2 { ' 

for all v in [0, 1]. It will then follow that 

X(0) = lim x(v) = ~ C B(0, v)dv > 0. (28) 

v->0 Jo 

This condition ensures that central shell-focusing singularity will be naked. 

We, now, discuss the existence of b(r) as a solution of the differential equation (26). We 
put 

x(r,v) = vb(r) + vh(r, v) + M(r,v) (29) 
a continuous function of r, in equation (26) and rearranging it, we get 

^ = ^{2B(r,v)x^]=f(x,r) (30) 
dr y/v 1 J 
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with the initial condition 



x(0, v) = vb(0) + vh(0, v) + M(0, v). (31) 

Let us ensure the existence of C l - function x(r, v) as a solution of above initial value 
problem defined throughout the cloud. The function f(x,r) is continuous in r, with x 
restricted to a bounded domain. With such domain of r and x, f(x, r) is also C 1 - function 
in x which means f(x,r) is Lipschitz continuous in x. Therefore, the differential equation 
(30) has a unique solution satisfying initial condition (31). 

Further, we can ensure that the solution will be defined over the entire cloud i.e. for 
all r in [0,r c ] by using the freedom in the choice of arbitrary function B(r,v). For this 
, we consider the domain [0, r c ] x [0, d] for some finite d. Let S be the supremum of the 
modulus of f(x,r). Then the differential equation (30) has a unique solution defined over 
the entire cloud provided [10] 

r c <inf(r c ,|) = |. (32) 



This yields 



max 

0<r<r c ,0<x<d 



' 2B(r,v)x 3/2 



< -• (33) 



Condition (32), in turn, will be satisfied if the weaker condition 

0<\B(r,v)\x 3 / 2 <^ (34) 

holds for all r in [0,r c ]. 

The collapsing cloud may start with r c small enough so that the expression d^/v/2r c 
which is always positive, satisfies the condition (34) with x restricted to a bounded domain. 
We then have infinitely many choices for the function B(r,v), which is continuous and 
satisfies conditions (27) and (34) for each choice of v. For each such B(r,v), there will be 
a unique solution x(r,v) of the differential equation (30), satisfying initial condition (31), 
defined over the entire cloud and in turn, there exist a unique function b{r) for each such 
choice of B(r,v), that is given by the expression 

b ^ = xjryu) - M(r,v)-vh(r,v) ^ 

over [0, r c ]. Thus, for a given constant v -surface and given initial data of mass function 
F(t,r) = r 3 M(r,v) and h(r,v) = (e r 92 ^ r ^ — l)/r 2 satisfying physically reasonable condi- 
tions (expressed on M), there exists infinitely many choices for the function b(r) such that 
condition (27) is satisfied. The condition continues to hold as v — > 0, because of continuity. 
Hence, the central singularity developed in the collapse is a naked singularity. 
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4 Definition of stability and genericity of NS/BH 



In this section, we give the appropriate definitions of stability and genericity. Our defini- 
tions of stability and genericity are based on the structural stability of a dynamical system, 
and genericity of a property of a dynamical system respectively. For more details, refer 
chapter 7, §3,4 of Abraham and Marsden [11]. 

Let M be a manifold on which a vector field or a dynamical system is defined. Let 
Af(M) be the space of all vector fields on M. Whitney C r -topology on X(M), generated 
by the norm is given by 



where V and W are vector spaces, U an open subset of V, and /:[/—> W. D k f denotes 
kth Frechet derivative of /. X(M) endowed with Whiteny C r - topology is denoted by 



Let Xi be a vector field or a dynamical system on M. Then X 1 is structurally stable if 
there is a neighbourhood $ of X\ G X r (M) in the Whitney C r - topology such that Y 1 G $ 
implies X\ and Y\ are topologically conjugate i.e. they have equivalent phase portraits. 
This means that there is a homeomorphism h Q : M — > M carrying oriented orbits of Xi to 
oriented orbits of Y\. 

A property of vector fields in X r (Wl) is a proposition P(xi) with a variable X\ G 
X r (M). A property P(x 1 ) with a variable X\ G X r {M) is generic if the subset {x\ G 
X r (M)/P(xi)} C X r (M) contains a residual set. 

A subset A of a topological space X\ is called residual if and only if A is the intersection 
of a countable family of open dense subsets of X\. A topological space Xi is a Baire space 
if and only if every residual set is dense. We also know that every complete metric space 
and in particular every Banach space is a Baire space. Also, whether M is compact or not, 
X r (M) is a Baire space. 

In our case, we apply these definitions analogously to the evolution of initial data into 
a gravitational collapse leading to a naked singularity or a black hole. We treat evolving 
initial data as a vector field or a dynamical system and consider the space of all initial data 
with sufficient differentiability defined on a collapsing compact spherical shell of Type I 
field, in place of X(M), endowed with a suitable C r -topology. Property P of a dynamical 
system becomes the property of initial data, namely, whether this initial data leads the 
collapse to a naked singularity or a black hole. Thus, the definitions of stability of a naked 
singularity and genericity of its occurrence can be stated as follows. 

Let I a be the initial data set, which when evolves, leads the collapse to a NS/BH. We 
say that a NS/BH is stable, if there is a neighbourhood I of I a in C r -topology such that, if 
Ii is another initial data in I, then I\ also leads the collapse to a NS /BH. In other words, if 
the set of initial data leading the collapse to a NS/BH forms an open subset of C r -space of 
all initial data, then the NS/BH will be stable. Similarly, occurrence of a NS/BH will be 
said to be generic, if the set of all initial data leading to NS/BH, is a dense subset of the 
parent C r -space. All sets of initial data are required to satisfy the constraint equations, in 
addition to the energy conditions. 




x r (M). 
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4.1 Stability of Naked Singularity 

The analysis in section III shows that the only conditions on the initial data which evolves 
the collapse into a naked singularity are the energy conditions. Hence, for stability, we 
have to examine energy conditions only What we now show is that, the set of such initial 
data satisfy energy conditions forms an open subset of the space of all initial data. For this 
, we use the technique of Saraykar and Ghate [12]. The energy conditions p>0; p + p r >0; 
P + Pi9>0 can be written in the form 

[3M + rM tr + rv'M )V ] = E l > 0, [3M + rM r - vM tV ] = E 2 > 0, 
[3wM+ rwM r + wiM tV - rv(v + rv')M )rv ] = E 3 > 

where w = 2(v +rv ')+rv v tT and w± = 2r 2 v' 2 — (2+r)v 2 v tT on the domain T>. Here, we have 
used an assumption that v is an increasing function of r so that v ' > on any surface t = tj 
say, so that no shell-crossing singularity condition (ie.R f > 0) holds. So, this is a valid 
assumption. We know that as the nuclear fuel exhausts in a star, it starts collapsing under 
its tremendous gravitational pull, and in this case, in the formation of central singularity 
radial pressure diverges, so we assume radial pressure p r to be non-negative throughout the 
gravitational collapse. Therefore, from equation (12), we have M„ < 0. Then, E\ < E 2 . 
Hence, only E\ and E$ shall take part in further stability analysis. 

We assume that X be an infinite dimensional Banach space of all C 2 real-valued func- 
tions defined on V, endowed with the norms 

|| M(r,v) || i= sup |M| + sup |M r | + sup |M„| 

and || M(r,v) \\ 2 = sup |M| + sup |M rr | + sup |M rv | + sup |M„„| 
v v v v 

These norms are equivalent to the standard C 1 and C 2 norms 

|| M(r,v) || 1= sup(|M| + |M r | + |M„|) 
v 

and || M(r,v) || 2 = sup(|Af] + |M rr | + | M rv | + |M W |) (36) 

v 

Let g = {M(r, v) : M > 0, M is C 2 , E x > and E 3 > on V} be a subset of X. 

We show that Q is an open subset of X. For Min Q, let us put S = min(M), 7 = min(£'i) 
, P = min(_E 3 ), A = max(w), Ai = max(tui), A2 = max(t> + rv') and A3 = max(t>') for 
varying r in [0, r c ] and ve[0, 1], the functions involved herein are all continuous functions 
defined on a compact domain V and hence, their maxima and minima exist. We define a 
positive real number 

= I • xa 2 JL JL_ JL JL JL i 1 

li 2 mmto, g , ^, 3 ^ Aa> m , 4 ^ A , 4Ai? 4rcA2 i- 

Let M x {r,v) be C 2 in D with || M— M x ||i< fi and || M— Mi || 2 < fi. Using definition 
(36), we get \M X - M| < // , |M 1;r - M r | < /i , |M M - M„| < // and |M 1;TO - M rt) | < // 
over V. Therefore, for choice of //, the respective inequalities are 

Mi > M- ^ > 0, 3|Mi -M\<1 
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r\M 1:r - M r | < r c \M 1>r - M r | < j, rv'\M hv - M„| < r c X 3 \M hv - M jV \ < 

3w\M 1 - M\ < f , rw\M 1>r - M r | < r c A|M 1>r - M r | < f , 

o o 

Wi|M M - M„| < r c Ai|Mi )t , - M„| < f , ru(u + n/)|Mi >rt) - M TO | < f (37) 

o o 

that are satisfied on T>. The II nd , III rd and 7T/* /l inequations from above yield 
3|Mi - M] + r|Mi >r - M r | + rV|M M - M >v \ < ^ < 7 < E 1 

Further, we can write | [3 Mi + rM^ r + rv'M ljV ] - E x \ < E 1 where E 1 > on £>. Hence, 
[3Mi + rM l r + rt> 'M 1; „] > on T>. Using similar analysis for last four inequations of 
equation (37), we obtain [3wMi +rwM^ r + wiMi >v — rv(v + rv')M 1>rv \ > on V provided 
w > and Wi > on T>. 

Thus, Mi > 0, Mi is C 2 , [3M 1 + rM 1>r + rv'M 1;V ] > and [3wMi + rwM l , r + WiMi^ - 
rv(v + rv')M ltrv ] > 011V provided w > and u>i > throughout V. Therefore, Mi(r, v) 
also lies in Q and hence, Q is an open subset of X. 

The conditions w > and u>i > required for openness of Q on P, impose conditions 
on the function u(t,r), v(t,r) = A{t,R) and A — J A^dt = r 2 g 2 (r,v), therefore, h(r,v) = 
e r 92{r,v)-i j r 2 j g cons trained by the conditions w > and w\ > in domain V. Now, let 
Y be an infinite dimensional Banach space of all C 1 real-valued functions defined on V. 
The set 7i of functions h(r, v) with restrictions w > and wi > would form (by similar 
technique) an open subset of Y. Then Q x 7i is an open subset in the product space X xY. 

Taking (Mi, hi(r, v)) in the neighbourhood of (M, h(r, v)) in G x H, and using equation 
(34) analogously for Mi and hi, we have a choice of infinitely many Bi(r,v), such that for 
each such -£>i(r, v), there will exist a unique &i(r) so that the initial data of mass function 
r 3 Mi and hi(r, v) together will lead the collapse to formation of a naked singularity. Thus, 
naked singularity arising from (M, h(r,v)) is C 2 - stable in the sense defined above. 

The analysis given above, guarantee the existence of a metric function v(t, r) for a given 
initial data set. Such choice of u(t,r) and expressions for G and H together will yield the 
metric (1) as an exact solution leading to the occurrence of naked singularity. Thus, given 
initial data and the weak energy conditions above determine the matter distribution and 
the metric of the space-time. 

4.2 Genericity of NS/BH 

Let us discuss the genericity of the occurrence of these singularities. We have seen that in 
Type I matter field collapse, the end state of collapse is governed by the choice of initial 
data. As explained above, occurrence of a NS is generic if the set Q x 7i of 'true' initial 
data ( M, h(r,v) ), contains a residual set in the space. The space X xY contains negative 
functions which cannot be realized as a limit of a sequence of positive functions belonging 
to Q x 7i. Hence, Q x TC is not dense in X x Y. The space under consideration is a Banach 
space. Every complete metric space and in particular every Banach space is a Baire space, 
and since in a Baire space, every residual set is dense, it follows that Q x 7i is not residual. 
Therefore, we conclude that occurrence of naked singularities is not C 2 -generic. Further, 
since the set Q x H is open, it cannot be nowhere dense and hence Q x Ti is not a meagre 
set either. Thus, it is a substantially "big" set. 
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5 Special case p r = p r (r) 

The case p r — 1 has been studied by Goswami and Joshi and it is shown that the presence 
of a non- vanishing pressure gradient gives rise to either formation of NS or a BH [13]. 
We illustrate above analysis by studying a special case where p r be a function of r alone. 
It gives a more clear picture of the occurrence of naked singularity and its stability. We 
choose two allowed free functions p r and u(t, r) as follows: 

Pr=p r (r), u(t,r)=c(t)+ri(R) (38) 

Now, let us see, how the Einstein's field equations react to this choice which will decide 
the evolution of collapse. As p r is a function r alone, integrating equation (4), we obtain 

F(t,r) = -^-R 3 + z(r) (39) 

where z(r) is another arbitrary function of r. We need to choose z(r), z(r) = |r 3 p r (r) so 
that the mass function takes the form 

F = ^(2r 3 -i? 3 ) (40) 
3 

that satisfies the regularity conditions at the initial epoch F(ti, r) = p r r 3 /3 and F(ti, 0) = 
0. Therefore, F > provided p r > 0. Hence, non- negativity of radial pressure shall 
maintain non- negativity of mass function throughout the gravitational collapse. Next, 
equation (3) becomes 

_ 1 3p r (2r 2 - R 2 R') + p' r (2r 3 - R 3 ) 
P ~ 3 R 2 R> { ' 

The density blows up at the singularity in the limit of approach to the singularity ie. as 
r — > and t — > t s . At the initial epoch t = U, 

P = p r + \rp' r > (42) 

, an all important relation between energy density and radial pressure at the initial epoch. 
Equation (6 ) takes the form 

G = e 2 ^ R) d p {r) (43) 

with d p (r) = 1 + r 2 b p (r) where b p is the energy distribution function for the collapsing 
shells. Using equation (38) in equation (5), we obtain 

(P + Pr) V,R R'R = 2(Pe - Pr)R' -Rp'r (44) 

Since the density blows up at the singularity, from above equation, it is clear that the 
tangential pressure also blows up in the limit of approach to the singularity. Further, we 
find that all of initial data is not independent at the initial epoch t = ti 

(p + Pr) V,R r = 2 (Pe - Pr) ~rp' r . (45) 

Above equation and equation (42) suggest that b p , radial and tangential pressures will be 
sufficient to form initial data set for the collapse. We are going to show the existence of b p 
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through p r and pg. Therefore, we consider p r ,pg as an initial data set. Using this initial 
data prescribed at the initial epoch t — ti to evolve the collapse, we can integrate the 
equation (44) and obtain, 

r,{R)=r 2{PB '7 Pr ' ) ~ RP ' r dR (46) 

/V ' Jo R(p +Pr o ) V ' 

Thus, velocity distribution function r)(R) is determined by smooth functions p r and pg. 
Using equations (40) and (43), equation (7) can be written as 



y/RR = - a{t) e v 01 + r 2 b p ]Re 2 ^ - R + j (2r 3 - R 3 ) (47) 

where a(t) is a function of time and we choose a(t) = 1 by suitable scaling of time co- 
ordinate. 
We define 

P 2»j(fl) _ 1 

h p( R ) = R 2 = 2 9p(R) + 0(R 2 ) (48) 
Using this definition, equation (47) takes the form 



y/v v = -sjvb p e^ + e 2 ^[v 3 h p {rv) + j (2 - v 3 )} (49) 
Integrating above equation, we obtain 

t(v,r) = / , V 50 

Jv ^Jve 4 % + e 2 ^[v 3 h p (rv) + f (2 - v 3 )} 

where the variable r is treated as a constant. Further, the time taken for the central shell 
to reach the singularity is given by 

= t — - (51) 



vb p (0)+v 3 h p (0) + f(0)(2 



v 



3\ 



and t So is well defined provided [ vb p (0) +v 3 h p (0) + p r (0)(2 - v 3 )/3 ] > 0. The time taken 
for other shells to reach the singularity can be determined from the Taylor expansion of 
t(v, r) around r = as t — > t s , 

t s (r) = t So +r Xp (v) + 0(r 2 ) (52) 

where Xp( v ) — dt(v, r)/dr\ r=0 and 

= _ 1 , ^(0) + ^V(0) + «!f(2-^] 
PV ; 2i„ [^ p ( )+t;3/ lp (o) + ^M(2-t;3)]3/2 v ; 

For any constant v surface, we have d(v) = which implies -Jv v — — (l/x P (v))\/v v' . 
Now, equation (49) takes the form 



yfiv' = X P (v)\ vb p (0) + v 3 h(0) + ^M( 2 - v 3 ) + 0(r 2 ) (54) 
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The radial null geodesic equation (24) in the variables (u = r a ,R) as t — > t s , r — > 0, gives 

*o /2 = ^X P (0)^ r (0) (55) 

From equation (53), it is clear that Xp( v ) depends on the functions b p (0), h p (0),p r (0) and 
these functions, in turn, depend on the given set of density, pressure and velocity function v. 
Therefore, equation (55) depicts that tangent to the singularity curve t = t So is completely 
determined by the given set of initial data and radial pressure in particular. The radial 
null geodesic emerging from the singularity in (R, u) co-ordinates is R = x Q u, while in (t, r) 
plane, the null geodesic near the singularity is given by t — t s (0) = x Q r 5 ^ 3 . 

It is clear that if x P (0) > then x Q > and with this, we get radially outgoing 
null geodesies coming out from the singularity, giving rise to a naked central singularity. 
However, no such trajectory comes out if x P (0) < as we have a black hole situation. 
Thus, we have some region of space-time for which x P (0) > depending on initial data, 
leading to occurrence of central naked singularity in that space-time region. Therefore, we 
choose bp to satisfy the following differential equation on a constant f-surface 

l ^[vy p (r)+v%, r (rv) + ^(2-v 3 )] J 

r^-^ dv — B(r,v) 56 

2 [vb p (r) + v 3 h p {rv) + ^(2 - v 3 )] 3 / 2 

where B(r, v) is a continuous function defined on V such that 

i >K(o) + A, r (o) + f(2-^)] , mn , n , 77 rnil ,, 7 , 

jr-r^- dv — B(0,v <0 for all v e 0, 1 . 57 

2 [vb p (0)+v 3 h p (0) + ^(2-v 3 )] 3 / 2 1 

It then follows that 

Xp (0) = limXpM = -f*B(0,v) dv > 0. (58) 

The above requirement ensures the nakedness of central shell focusing singularity. Now, 
we discuss existence of b p as a solution of the differential equation (56). For this, we use 

y(r, v) = v b p {r) + v 3 h p {rv) + ^1(2 - v 3 ) (59) 
a continuous function of r, in equation (56), and we obtain 

^. = ±[2B(r-v)y 3 / 2 } =<p{y,r) (60) 



with initial condition 



y(0, v) = vb p (0) + v%(0) + PM(2-v 3 ). (61) 

This initial value problem defined throughout the cloud has a C 1 - function y(r,v) as its 
unique solution defined over the entire cloud. Using the arguments of general case (section 
III) for the function <f(y,r) and B(r,v) ) a weaker condition 

0<\B(r,v)\y 3 ^<^ (62) 
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holds for all r in [0, r c ] and v G [0, 1] where n is a finite number such that < y < n. 
This condition holds true for y restricted to a bounded domain and the collapse may start 
with r c small enough. Therefore, there are infinitely many choices for the function B(r,v), 
which is continuous and satisfies the conditions (57) and (62) for each choice of v. For 
each such £>(r, v), there will be a unique solution y(r,v) of the differential equation (60), 
satisfying initial condition (61), defined over the entire cloud. This unique solution y(r,v) 
for each such choice of B(r,v) yields 

y(r,v) -v 3 hJrv) - ^(2-v 3 ) . . 

b p = ^-J. ^— j (63) 

over [0, r c ]. As each function appearing in the above expression is continuous, therefore, the 
condition holds true as v — > 0, because of continuity. We know that velocity distribution 
function i](R) is determined through smooth functions p r , pg on a given constant -y-surface, 
which in turn, through equation (48) determines the function h p (rv). Thus, for a given 
constant v -surface and given initial data of radial and tangential pressures p r and pg, there 
exists infinitely many choices of the function b p such that condition (57) is satisfied. We 
have thus proved the existence of energy distribution function b p for a given set of initial 
data leading to the occurrence of a naked central shell focusing singularity. 

Now, we illustrate the stability of occurrence of NS. For simplicity, we consider the 
evolution of initial data on the initial surface t — ti where v — 1 and v' — 0. Using 
equations (42) and (45), the energy conditions p>0; p + p r >0; p + pg>0 can be described 
as 

P = Pr + \rp' r = K x > 0, 2p r + Kp' r = K 2 > 0, 

2(2 + rr] tR )p r + [1 + 1(2 + r^ R )]rp' r = K 3 > 

We denote by X p be an infinite dimensional Banach space of all C 1 real-valued functions 
defined on [0, r c ], endowed with the C 1 norm. Let us take G p = {p T (j') '■ Pr > 0,p r is 
C 1 ,K 1 > and K 3 > on [0, r c ] } be a subset of X p . 

Now, to show that G p is an open set in X p . For p r in G p , we put a\ = mm(p r ), 
a 2 = min(if 1 ), a 3 = mm(K 3 ), and \ p = max(2 + r rj^). We define a positive real number 

^ _ 1 m j n | Q , a 2 a 2 «3 «3 «3 | 

2 3 t c 6Ap 3?^ c 3v c \ p 

Let p ri (r) be C 1 in [0, r c ] with || p r — p ri \\< p. Now, using definition of standard C 1 
norm, we get \p r — p ri \ < p and \p' r — p' \ < p over [0, r c ]. Then, for each choice of p, the 
respective inequalities 

OL\ | . a 2 1 , 1 , / i «2 

Pri>Pr-y>0, |p n -PH<y, -Pr\ < g^bn "Pri < "g-> 

2(2 + r^ R )|p ri -p r | < 2A p |p ri -p r | < y, r|p r / < r c |p^ < y, 

T Ct 3 

-(2 + rri, R )\p' ri -p r \ <r c X p \p' ri - p' r \ < — 



14 



are satisfied in [0,r c ]. Using above inequalities analogous to the general case, we obtain 
p ri > 0, [p ri + rp'jS ] > 0, [2 (2 + rr] :R )p ri + [1 + (2 + rr}, R )/3 }rp' ri ] > 0, provided 
(2 + ri]^) > on [0, r c ]. Therefore, p ri G G p and hence, G p is an open set in X p . 

The function r)(R) is determined by p r , po on the given constant v -surface and the 
condition (2 + rr] >R ) > is necessary for openness of the set G p . Therefore, the C 1 function 
h p (R) = (e 2T,( ^ — l)/R 2 is also determined by initial data p r , p e on a given constant 
v -surface constrained with (2 + rr] tR ) > on [0, r c ]. Let 7i p be a subspace of X p containing 
such initial data h p (R). Then G p x 7i p is an open set in the product space X p x X p . 

Let (p ri , h pi (r)) be in the neighbourhood of (p r , h p (r)) in G p x H p . Using equation (62) 
analogously for p ri and h pi (r), we have a choice of infinitely many Bi(r, 1), such that for 
each such Bi, there will exist a unique &i(r) so that the initial data p ri and h pi (r) together 
will lead the collapse to formation of a naked singularity. Thus, naked singularity arising 
from initial data (p r ,h p (r)) is C 1 - stable. Thus, this particular case illustrates clearly 
how the choice of initial data leads the collapse to a naked singularity and how it is stable 
under small perturbations of initial data in an appropriate mathematical sense. 

6 Discussion and conclusions 

1. Our main conclusions in the paper are the following: (a) Given a C 2 - mass function 
M(r,v) and a C 1 - function h(r,v), on any v = constant surface, we find a C 1 - energy 
distribution function b(r) such that the collapse ends in a naked singularity, (b) With 
physically reasonable conditions put on mass function M, and with C 1 - function h(r,v), 
the initial data consisting of (M, h(r,v)) leading the collapse to a naked singularity forms 
an open subset of the space X x Y . This establishes stability of naked singularity with 
respect to initial data. 

2. Similar analysis for the case x(0) < ( that is B(0,v) > for all v in [0, 1]) shows 
that the set G x H of 'true' initial data (M, h(r,v)), leading the collapse to a black hole 
forms an open set and it is neither dense nor nowhere dense in the space of all initial data. 
Thus, occurrence of black hole is stable but not generic. 

3. Our analysis is carried out with given initial data on a constant v- surface. Hence, 
the role of Mand b(r) can be interchanged (refer equation 35) such that the mass function 
would evolve in the continuous gravitational collapse. 

4. The conditions w > and w 1 > which occurred in our stability analysis of NS 
(general case) can yield some information about the function v. For simplicity, let us 
consider these conditions on the initial surface t = ti. For w > 0, we obtain z/ r > — 2/r for 
all r € [0, r c ]. This means z/ r G (— 2/r c , oo). Next, for w\ > 0, we get [—(2 +r)z/ r ] > 0, 
this implies v r e (— oo,0). Thus, common domain of v T is (— 2/r c ,0). Therefore, v is a 
decreasing function of r in [0,r c ] for sufficiently small r c . 

While, in particular case [2 + rr) >R ] > yields r\ r e (— 2/r c , oo), therefore, rj(r) is 
a decreasing function of r in [0, r c ]. Also, the above said condition can be written as 
[Q(p r + Pe) — fp'r] > 0, it means C 1 - stability of occurrence of NS is associated with initial 
data p r and p e such that [6(p r + pe) — rp' r ] > 0. It is this choice of data that gives rise to 
stable NS. 

5. The space X used in the analysis is an infinite-dimensional Banach space, and as 
such, is not locally compact. It is not therefore, straightforward to define a measure on 
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this space. Hence, we can not say whether the set Q is of zero measure or not. However, as 
mentioned earlier, Q is not meagre, and thus, it is a substantially 'big' set. Similar argument 
is applicable to the set Q x 7i in the product space X xY . We believe that rigorous analysis 
that we have worked out shall help further in understanding the gravitational collapse of 
Type I matter fields. 
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